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Abstract. This note considers cooperative irreducible systems of ordinary 
differential equations admitting a C 1 first integral with positive gradient. We 
prove that all forward or backward nonwandering points are equilibria. We 
obtain also some results on the global phase portrait of such systems. The 
main tool of proof is a canonically denned Finsler structure with respect to 
which the derivative skew-product dynamical system is contractive. 



1. Introduction 
A system of ordinary differential equations (ODEs) 
(I) ± l = f(x), x = (x\...,x n ), 

where / = (f 1 , ...,/"): X -> K" is a C 1 vector field on an open set X C R", is 
called cooperative if (df 1 /dx 3 ) (x) > for ij^j and all x G X. 

Let <p(t;xo) denote the nonextendible solution of system (TTJ) with the initial 
condition 0(0; xq) = xq. We write (f>tXo instead of </>(_; xo). For each x G X 
the mapping t i— > cf> t x (called the trajectory of x) is defined on an open interval 
(o(x), t{x)) containing 0. The restriction of the trajectory of x to (c(x), 0] [resp. to 
[0,t(x))] is called the backward [resp. forward] semitrajectory of x. The images of 
(semi)trajectories are referred to as (semi) orbits . We say x G X is an equilibrium if 
4>tx = x for all tgl, or, equivalently, if f(x) = 0. A point y G X is an u-limit point 
of x G X if there is a sequence t/. — > oo as k — > oo, such that limfc._i.oo 4>t k x = y. 
Notice that r(x) = oo, while it is possible that r(y) < oo. The definition of an 
a-limit point is analogous. The set of w-limit [resp. a-limit] points is called the 
u) -limit set [a-limit set] of x, and denoted by uj(x) [a(x)\. 

The symbol ||-|| stands for the Euclidean norm in K™. We say x G X is forward 
nonwandering if for each e > and each < t < t(x) there are y G X with r(y) > t 
and t < 9 < r(y) such that ||x — y\\ < e and ||x — (t>ey\\ < e. A point x G X is 
backward nonwandering if for each e > and each o~(x) < t < there are y G X 
with <r(y) < t, and o~(y) < 9 < t such that ||x — y\\ < e and ||x — <frey\\ < e. A point 
that is forward nonwandering or backward nonwandering is called nonwandering . 
It is straightforward that an w-limit point is forward nonwandering. Notice that in 
the above definitions we do not assume the forward [backward] semitrajectory of 
either x or y to be defined on the whole half-line [0, oo) [(— oo, 0]]. 
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For two points x, j/£R" denote 

x < y if x % < y l for each i, 
x < y if x < y and x ^ y, 
a; < y if x % < y l for each i. 

For x < y we define a closed order interval as 

[as,y] :={zeR n :x<z<y}, 

and for x <C y we define an open order interval as 

[[ac, 2/]] :={zel":i«z« y}. 

A set X C 1" is said to be p-convex if the line segment with endpoints x and y is 
contained in X for each x, y G X, x < y, and order-convex if [a;, j/] C X for each x, 
y G X, a; < y. 

The next result gives an important property of cooperative systems of ODEs. 

Theorem 1. Assume ^ is a cooperative system of ODEs on a p-convex open set 
X C M". Let x < y. Then fax < fay for each t G [0, min(r(a;) ) r(j/))). 



Miiller (1927)] an d [Kamke (1932)] . Some 
Wazewski (1950)] . The property is referred 



The above theorem was proved in 
gaps in the earlier proofs were filled in 
to as monotonicity of the local flow generated by ([1]) 

An important feature of cooperative systems of n ordinary differential equations 
is that the limiting behavior of a point whose forward semiorbit has compact closure 
is at most so complicated as that in a general system of n— 1 ODEs. More precisely, 
the following result holds (see Theorem A in Hirsch (1982)] ) : 

Theorem 2. Let ^ be a cooperative system of ODEs on a p-convex open set 
X C M". Assume that the forward [resp. backward] semiorbit of x G X has compact 
closure in X. Put L = lj(x) [resp. L = a(x)]. Then the restricted flow {fa\L} 
is topologically equivalent to the flow of a Lipschitz system of ODEs on 
Moreover, no two points in L are related by <C- 

An important class of cooperative systems is formed by cooperative irreducible 
systems of ODEs, that is, cooperative systems such that for each x G X the matrix 
[(df l / dx 3 )(x)\ is irreducible. A system of ODEs is called strongly cooperative if 
(df/dx j )(x) > for i ^ j. 

For cooperative irreducible systems Theorem [T] can be strengthened to: 

Theorem 3. Assume |lj is a cooperative irreducible system of ODEs on a p-convex 
open set X C K™. Let x < y. Then <j) t x ^ fay for each t G (0, min(r(x), r(y))). 

The above property is called strong monotonicity of the local flow generated by 
©■ 

For cooperative irreducible systems one can say much more about their behavior 
(see Theorem 2.4 in |Smith and Thieme (1991*)] ; for earlier results see [Hirsch (1988)] 
and |Polacik (1989)] ): 

Theorem 4. Assume (JlJ is a cooperative irreducible system of ODEs on a p-convex 
open set X C W l such that each forward semiorbit has compact closure in X . Then 
there exists an open dense set Y G X such that oj(x) is a singleton for each x G Y. 
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Generally, in cooperative irreducible systems there are no restrictions on a-limit 
sets. More precisely, according to Smale (1976)] , any dynamics on the standard 



n-dimensional simplex can be embedded as a repeller in an (n + l)-dimensional 
strongly cooperative system of ODEs. 

For a recent monograph on cooperative systems of ODEs sec [Smith (1995)] . 



2. Cooperative Irreducible Systems with First Integral 

By a first integral for ([1]) we mean a continuous function H : X — > K which is 
constant on orbits of ([T]). A first integral is nontrivial if it is not constant on any 
open set. 

The existence of nontrivial first integrals puts severe restrictions on cooperative 
irreducible systems. Namely, it was proved in Hirsch (1985)) that if the set of 



equilibria is countable and all points have forward semiorbit closure compact in X 
then each first integral is trivial. 

In the three-dimensional case much more can be proved, even without assuming 
the abundance of points with compact forward semiorbit closure (see 
Mierczyhski (1995)] ): 



Theorem 5. Assume that (TTJ) is a cooperative irreducible system of ODEs on R™ 
admitting a C 1 first integral with nonzero gradient. Then each limit set is either 
empty or a singleton. 



3. Cooperative Irreducible Systems with Monotone First Integral 

3.1. Limiting Behavior. When one assumes that a first integral for a cooperative 
irreducible system of ODEs is strongly monotone, that is, from x < y it follows 
H{x) < H(y), then, under the assumption that all forward semiorbits have compact 
closure in X, all w-limit sets are singletons (compare |Mierczyhski (1987)] ). 

This result carries over to the case of abstract strongly monotone semidynamical 
systems with monotone first integral defined on strongly ordered Banach spaces 
( [Arino (1991)] ). 

Further, for some classes of cooperative periodic [resp. almost periodic] (in 
time) systems of ODEs admitting a first integral with appropriate monotonicity 
properties it was proved that each solution with compact forward semiorbit clo- 
sure converges to a periodic [resp. almost periodic] solution ( [Nakajima (1979)] , 
Sell and Nakajima (1980)] , [Tang et al (1993)] , [Jiang (1995)] ). Moreover, in the 
almost periodic case in the corresponding nonlinear skew-product (local) flow the 
image of that solution intersects each fiber at precisely one point (in the case where 
the system is periodic with period T this simply means that the limiting solution 
has period T). 

In many of the proofs of the results mentioned above the idea was to use a kind 
of Lyapunov function. 

In the present note we do not assume any compactness of forward (or backward) 
semiorbits. On the other hand, we make extensive use of the fact that a cooperative 
irreducible system of ODEs generates a linear skew-product dynamical system on 
the tangent bundle of X possessing some monotonicity properties. 

Let us introduce some notation. For t G (cr(x),r(x)) the derivative D<fi t (x) of 
4>t{x) with respect to X is a linear isomorphism from the tangent space at x into 



4 



JANUSZ MIERCZYNSKI 



the tangent space at <f>tx, satisfying the nonautonomous linear matrix ODE 

M' = Df{cj) t x)M 

with initial condition M(0) = Id, where Df := [(dp/dx^)]. The local linear skew- 
product dynamical system 

(x, v) H> ((f> t x, D<j) t (x)v), x & X, v e M™, a(x) <t< t(x), 

will be referred to as the derivative local flow. 

The order relations <, < and <C are defined in a natural way on tangent vectors. 
We will refer to vectors v > as nonnegative, and to vectors u > as positive. 
The set of all nonnegative (free) vectors is called the (nonnegative) cone. 

The derivative flow enjoys the following strong monotonicity property (see 
Hirsch (1985)] ): 



Theorem 6. Assume that §Q is a cooperative irreducible system of ODEs on an 
open set X C M. n . Then for each x G X . each t G (0, t(x)) and each v > one has 
Dcj) t (x)v > 0. 

Let H(x) stand for the level set of the first integral H passing through x. If H 
is of class C 1 and grad-ff(x) ^ for all x 6 X then H(x) is a C 1 submanifold of 
codimension 1. 

The main result of the present note is 

Theorem 7. Let §Q be a cooperative irreducible system of ODEs on an open 
X C R n . Assume that ([1]) admits a first integral H of class C 1 with positive 
gradient. Then each nonwandering point is an equilibrium. 



Proof. In Mierczyhski (1991)] it was proved that for a cooperative irreducible sys- 



tem of ODEs admitting a C 1 first integral with positive gradient there is a canonical 
Finsler on the foliation of X into level sets of H under which the derivative (local) 
flow is contractive. By a Finsler we understand a continuous mapping (x, v) i— > \v\ x 
such that for each x the assignment v t— > \v\ x is a norm on the tangent space at x 
oiH(x). 

This canonical Finsler is constructed in the following way: For each x G X, 
we translate the tangent space of T-l(x) at i by a positive vector v x such that 
(grad-ff(x), v x) — 1< The intersection of the resulting hyperplane with the nonneg- 
ative cone is a compact convex set A x linearly isomorphic to the standard (n — 1)- 
dimensional simplex. Finally, we take the compact convex balanced set A x — A x to 
be the unit ball in the tangent space of H(x) at x. 

Proposition 2 in Mierczyhski (199T)] states that for each x G X, each nonzero 
vector v tangent at x to %(x) and each t G (0,t(x)) one has \D(j) t (x)v\^ t ^ < \v\ x . 
(As a matter of fact, that proposition is stated for strongly cooperative systems, 
but its proof carries over verbatim to the case of cooperative irreducible systems.) 

Suppose for contradiction that a; is a forward nonwandering point not being an 
equilibrium. Let L C X be a C 1 embedded (n — l)-dimensional disk transverse 
to f(x) and having x in its relative interior. Pick s > 0, s < min{r(z) : z G L}. 
Denote 

A := max{\D(l> s (z)v\ < f >s ( z ) : z G L, v G T Z H, \v\ z = 1}, 

where T Z TL denotes the tangent space of 'H(z) at z. It is obvious that < A < 1. 
Now, take a C 1 embedded (n — l)-dimensional disk M C L transverse to f(x), 
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having x in its relative interior and such that 

(1 - »)\f(x)\ x < \f(z)\ z < (1 + for all z G M, 

where /i := (1 - A)/(l + A). 

As x is, by assumption, a forward nonwandering point, there exists a point z G M 
and i > s such that </>tZ G M . But 

|/(<Mk* = |-D^t-.(0.*)/(MU* < l/(Mk« < < (1 - A»)|/(ar)U, 

a contradiction. 

For a backward nonwandering point x not being an equilibrium, let L C X be 
a C 1 embedded (n — l)-dimensional disk transverse to fix) and having x in its 
relative interior. Pick s < 0, s < max{cr(z) : z G L}. Denote 

A := mm.{\D^ s (z)v\^ z -) : z G L,v G T Z H, \v\ z = 1}. 

It is obvious that A > 1. Now, take a C 1 embedded (rt — l)-dimensional disk M C L 
transverse to f(x), having x in its relative interior and such that 

(1 - n)\f{x)\ x < \f(z)\ z < (1 + n)\f(x)\ a for all z e M, 

where fj, := (A - 1)/(A + 1). 

As x is, by assumption, a backward nonwandering point, there exists a point 
z G M and t < s such that (f>tz G M. But 

= |i?&- s (<M/(<MI<m > > A|/(*)|, > (i - m)I/(*)U, 

a contradiction. □ 



As a corollary we obtain (see Mierczyhski (1991) ) 



Theorem 8. Let (|T|) be a cooperative irreducible system of ODEs on an open 
X C K." admitting a first integral H of class C 1 with positive gradient. Then 

(i) Each u> -limit set is either a singleton or empty. 

(ii) If a(x) is nonempty then x is an equilibrium. 

Proof. Part (i) follows by the fact that each w-limit point is forward nonwander- 
ing. Also, each a-limit point is backward nonwandering. Assume that for some x, 
a ( x ) = {y}- We have = \f(y)\y = limt^-oolfj&xM ttx = sup{|/((/» t a;)|^ ta! : t G 
(— oo, t(x))} (see Proposition 2 in [Mierczyhski (1991)] ), hence |/(x)|a; = 0, and x 
is an equilibrium, therefore x = y. □ 

3.2. Global Picture. In the present subsection we shall obtain some insight into 
the global nature of the dynamical system restricted to a level set H(x) of H . The 
standing assumption will be: 

X C W 1 is an open order-convex set such that for any pair x, y G X their 
maximum x V y and minimum x Ay are in X . 

For instance, X = R™, or X = [[a, b]] for some a, b G R", a <C b, or else X is the 
positive orthant R™ := {x G R™ : x i > for all i}. 

We begin with an auxiliary 

Lemma 9. Let H be a C 1 first integral with positive gradient for a cooperative 
irreducible system {!]). Then for each x G X, H{x) is connected. 
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Proof. Fix x G X, and write Ti '■= Ti{x). Take y, z G Ti, y ^ z. As y + z = 
y Az+yV z, it follows that the points y, z, yAz and yVz belong to a two-dimensional 
affine subspace V. We have -ff(y A z) < H(x) and H(y V z) > H(x). For each 
A G [0, 1] denote by A> the union of the line segment joining y\f z with Xy + (1 — X)z 
and the line segment joining yAz with Xy + (1 — A)z. The set A := U Ae j l] ^ A 
is a two-dimensional (analytic) submanifold-with-corners contained in V. It is 
apparent that the gradient of the restriction H\A is everywhere nonzero. The 
implicit function theorem yields that the set A n Ti is C 1 diffeomorphic to the real 
interval [0,1]. This finishes the proof. □ 



As in the case of a Riemannian metric, for a C 1 curve 7 contained in Ti define 
the length of 7 as 

6 

a 

where 7 : [a, 6] — > % is a parametrization of the curve. It is straightforward that 
the length does not depend on a parametrization. We define the Finsler distance 
d{x, y) between two points x, y G Ti as the infimum of the lengths of all C 1 curves 
with endpoints x and y. The C 1 manifold Ti together with the Finsler distance 
■) is a metric space. 

Theorem 10. Let a cooperative irreducible system (JlJ admit a first integral H 
of class C 1 with positive gradient. Assume that T-L is a level set of H such that 
t(x) = 00 for all x G W. Then either 

(a) There is precisely one equilibrium y in %, and y is a global attractor in T~L; 

or 

(b) There is no equilibrium in H, and for each x G H one has w(x) = 0. 

Proof. As a consequence of Theorem[71 if there is no equilibrium in % then u)(x) = 
for all x G 7i. So, assume y G % is an equilibrium. First, we claim that y is a 
unique equilibrium in %. Suppose per contra that there is another equilibrium y\. 
Put z := y V y\. We have z > y, z > y%. By strong monotonicity, for each t, 
< t < t(z), one has (j> t z <p t y = y and faz ^> 4>tVi = Vi- But this implies that 
(j) t z 3> z, that is, in the level set of H passing through z there are two points being 
in the <C relation, which is impossible. 

By Main Theorem in Mierczyhski (1991)] , y is (locally) exponentially asymp- 



totically stable relative to %. Consequently, the set A := {x G H : oj(x) = {y}} 
is relatively open in TL. Suppose by way of contradiction that A ^ H. Since by 
Lemma [9] % is connected, the relative boundary bd% A is nonempty. Pick a point 
z G bd« A. As z ^ A, we have u)(z) — by Theorem[7] 

Because y is asymptotically stable in H, by Conley (1978)) there is a compact 



relative neighborhood B of y in T-L such that B C A and </> f B C £> for all t > 0. Take 
e > so small that {x G Ti : <i(a;, y) < 2e} is contained in the relative interior of 
B. Pick 4 with z) < e. Let T > be such that d(4>rxi, y) < e. As cjypz 
exists, we must have d{<f>T%i,4>Tz) < £■ Consequently, d(4>TZ,y) < 2e. But from 
this it follows that 4>tz G B C A, hence w(z) = {y}. This contradiction completes 
the proof. □ 
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If X = R™, a well-known condition guaranteeing r(x) = oo for each x is the 
existence of positive constants C\ and C2 such that ||/(x)|| < Ci||x|| + Ci. For 
another condition see the following result. 

Theorem 11. Let a cooperative irreducible system JT]) defined on K™ admit a first 
integral H of class C , such that all the coordinates o/gradiJ(x) are positive, 
bounded and bounded away from zero, uniformly in x £ X . Then we have either 

(a) For each x £ R™, lo{x) is a singleton. Moreover, the set of equilibria is 
simply ordered by <C; 

or 

(b) For each x £ M™, uj{x) — 0. 

Proof. We begin by showing that there is a constant C > such that ||u|| < C\v\ x 
for each x £ R n and each vector v tangent at x to TL(x). Fix x £ R™, and put 

gradff(x) 
{ ' ' ||gradF(x)|| 2 ' 

Recall that in the construction of the Finsler | • | the unit ball B x is defined as A x — A x , 
where A x :— {v > : (G(x),v) = 1}. In particular, G(x) £ A x . Fix v £ B x , that is, 
\v\ x < 1. Write w = t>i — V2, where «i, 1)2 6 ^l^. For i = 1, 2, put ruj := Vi — G(x). 
Of course, v — w± — w-i- As Vi > 0, we have — uii — G(x) — Vi < G(x). On 
the other hand, (G(x), —Wi) is easily seen to be zero. Write G(x) — (01, . . . , a n ), 
—wi — (61, . . . , b n ). We have aj > 0, bj < aj and 2j=i = 0- Dchnc 

1 

Cj : " a,||gradiJ(x)|| 2 ~ a ' ? '- 

It is straightforward that all Cj's are positive and bounded uniformly in x £ R". 
We claim that 

- Cj < bj < c 3 

for each 1 < j < n. Indeed, suppose b^ > Ck for some k. We have then 



> ||gradg(x)ll 2 ~ °' = l|G(x)l ' 2 ~ a ' = - H a o h ^ 



j=t 3=1 
3#fe j¥k 

hence X)j=i a j^i > 0j a contradiction. The other inequality is proved by a similar 
argument. We have thus obtained that < ||tOi|| + \\w2W does not exceed a 
constant independent of x £ R n . 

Now suppose by way of contradiction that T :— t(x) < 00 for some x £ R™. 
Then the improper integral Jq \f{4>tx)\ ( f >tX dt is convergent, from which it follows in 
a standard way that the finite limit lim t _).T 4>t% exists, a contradiction. 

Assume that there exists an equilibrium y £ R™. We may assume y — and 
H(0) = 0. Take a positive real number r. Let e > be such that for each x 
with ||x|| < e one has H(x) < r/2. Since the coefficients of gradTJ are positive and 
bounded away from zero, for each x > 0, ||x|| = e, the half line {sx : s > 1} intersects 
the level set 7f^ 1 (r) := {z £ R™ : H(z) — r} at precisely one point M(x). It is 
easy to see that the mapping M is a homeomorphism of {x £ R n : x > 0, ||x|| = e} 
onto {z £ i? _1 (r) : z > 0}. As the latter set is forward invariant, a well-known 
application of Brouwer's fixed point theorem implies that there exists an equilibrium 
in i/ _1 (r). An analogous argument applies to the case r < 0. From Theorem [TUl we 
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deduce that w(x) is a singleton for any x € M. n . The fact that the set of equilibria 
is simply ordered follows by Proposition 2.1 in Mierczyhski (1987)] . □ 
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